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& . ABSTRACT 

^ ' Wc investigate the stellar populations of Lyman a emitters (LAEs) at z = 3.1 and 

3.7 in 0.65 deg 2 of the Subaru/ XMM- Newton Deep Field, based on rest-frame UV- 
to-optical photometry obtained from the Sub&ru/ XMM-Newton Deep Survey, the 
UKIDSS/Ultra Deep Survey, and the Spitzer legacy survey of the UKIDSS/UDS. 
Among a total of 302 LAEs (224 for z = 3.1 and 78 for z = 3.7), only 11 are detected 
in the A band, i.e., brighter than A'(3er) = 24.1 mag. Eight of the 11 A-detected 
LAEs are spectroscopically confirmed. In our stellar population analysis, we treat K- 
■ detected objects individually, while A-undetected objects are stacked at each redshift. 

We find that the A-undetected objects, which should closely represent the LAE pop- 



ulation as a whole, have low stellar masses of ~ 10 s - 10 8 5 M M , modest SFRs of 1 
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100 M© yr -1 , and modest dust extinction of E(B - V)+ < 0.2. The A-detected ob- 
jects are massive, Af star ~ 10 9 - 1O 1O5 M0, and have significant dust extinction with a 
median of E(B — V)* ~ 0.3. Four A-detected objects with the reddest spectral energy 
distributions, two of which are spectroscopically confirmed, are heavily obscured with 
E(B — V)* ~ 0.65, and their continua resemble those of some local ULIRGs. Inter- 
estingly, they have large Lyman a equivalent widths ~ 70 - 250A. If these four are 
d ' excluded, our sample has a weak anti-correlation between Lya equivalent width and 

Afstar- We compare the stellar masses and the specific star formation rates (sSFR) of 
LAEs with those of Lyman-break galaxies (LBGs), distant red galaxies, submillimetre 
galaxies, and I- or A-selected galaxies with photometric redshifts of z p hot ~ 3. We 
find that the LAE population is the least massive among all the galaxy populations 
in question, but with relatively high sSFRs, while NIR-detected LAEs have M star and 
sSFR similar to LBGs. Our reddest four LAEs have very high sSFRs in spite of large 
Afstar, thus occupying a unique region in the M sta r versus sSFR space. 

Key words: cosmology: observations — galaxies: formation — galaxies: evolution 
— galaxies: high-rcdshift — galaxies: stellar content — 



1 INTRODUCTION 

Lyman a emitters (LAEs) are a galaxy population which are 
common in the high redshift Universe. In the last decade, 
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many observations have succeeded in detecting LAEs 
from z ~ 2 up to z ~ 7, primarily based on narrow-band 
imaging to isolate Lyman a emiss ion (e.g. iHu et al l 119981 : 
iRhoads et al]|2000l : live et alj |2006j). Over a thousand LAEs 



Hu et a 


.1 120021: lOuchi et all 120031: iMalhotra & Rhoadd 


20041; 


Taniguchi et all 120051; Shimasaku et al. 


120061; 


Kashikawa et al. 2006; Dawson ct al. 2007; Muravama et all 


2007; Gronwall ct al. 2007: Ouchi ct al. 20Ql). 



Most LAEs have relatively faint, blue UV continua and 
large Lya equivalent widths (EWs). These properties col- 
lectively suggest that they are young star-forming ga lax- 
ies with low metallicities (e.g. IMalhotra fe~R hoads 2002|). In 
this sense, they may serve as building blocks of larger galax- 
ies in hierarchical galaxy formation. 

Stellar population analysis is helpful to further con- 
strain the nature of LAEs. Recent multiwavelength observa- 
tions covering near-infrared wavelengths have enabled anal- 
ysis of the stellar populations of LAEs, and there is grow- 
ing evidence that not all LAEs are primordial as described 
above. Table [T] summarises studies of the stell ar populations 
of LA Es, including the work presented here. iGawiser et all 
(|2007l ) have performed a stacking analysis of 52 LAEs at 
z — 3.1, and concluded that LAEs have low stellar masses 
(~ 1O 9 M0), young-age components ( ~ 20 Myr), and small 
dust extinctions (Ay = 0) (see also Gawiser et all 120061 ; 
iNilsson et ai1l2007h . iPirzkal et ail (|2007t ) have studied nine 
LAEs at z ~ 5 found by HST/ACS slitless spectroscopy in 
the Hubble Ultra Deep Field, to show that these faint LAEs 
are all very young (~ 1 - 20 Myr) with low masses (~ 10 6 

- 10 8 Mq) and small dus t extin ctions (Ay = - 0.6). On 
the other hand. lliai et all |2007l ) have concentrated on three 
spectroscopically confirmed LAEs at 2 = 5.7 with IRAC de- 
tection, and derived relatively large stellar masses, ~ 10 9 

- lO 1(, A/ , mild dust extinctions (E(B - V) ~ 0.15 - 0.2) 
and a wide range of age, ~ 5 - 100 Myr. lLai et all (|2008T ) 
have studied both IRAC-detected and undetected LAEs at 
z = 3.1 by the stacking method, to find that LAEs posses 
wide ranges of age (~ 160 Myr - 1. 6 Gyr) and mass (~ 10 8 - 
1O 1O M ) without dust extinction. iFinkelstein et all (|2009h 
have analysed 14 LAEs individually, which are either IRAC- 
detected or undetected. Their results also show wide ranges 
of stellar population age (~ 2.5 - 500 Myr), stellar mass 
(~ 10 s - 10 10 M o ), and dust extinction (A1200 ^ 0.3 - 5.0). 

However, most of the samples constructed to date are 
not necessarily large and deep enough to constrain the aver- 
age stellar populations of LAEs and to study rare, massive 
LAEs which may be a bridge to more massive and/or evolved 
objects like LBGs. In addition, correlations between stellar 
population parameters have not been addressed well. 

Recently, lOuchi et all (|2008T l have constructed the 
largest available sample of 2 = 3.1 and 3.7 LAEs in an about 
1 deg 2 of the Subaru/ XMM-Newton Deep Field (SXDF) 
from deep optical broadband and narrowband data. These 
large survey data enable us not only to find many mas- 
sive LAEs, but also to place better constraints on less 
massive (i.e., average) LAEs using stacking analysis. Our 
sample used in this study consists of 224 (78) LAEs at 
z = 3.1 (3.7) covered by deep JHK images taken with 
the UKIRT / WFC AM fro m UKIDSS Ultra Deep Survey 
jUDS; Warren et all [2007b and 3.6 - 8.0/Ltm images taken 
with the Spitzer/IRAC from the Spitzer legacy survey of 



the UDS field (SpUDS; PI: J. Dunlopfl among which 5 
(6) are detected in the K band (i.e., brighter than the 3<r 
detection limit of the if -band image). The K band corre- 
sponds to the rest-frame V and B bands at z = 3.1 and 
3.7, respectively. For Tf-undetected LAEs we make median 
stacked images for individual bandpasses. Since the vast ma- 
jority (96%) are undetected in K, the spectral energy dis- 
tributions (SEDs) constructed from stacking should closely 
represent the whole LAE population at each redshift. Since 
the time interval between the two redshifts is not so large 
(the age of the universe is ~ 2.0 Gyr at z = 3.1 and ~ 1.7 
Gyr for 2 = 3.7), we treat the LAEs at 2 = 3.1 and 3.7 
collectively as objects at 2 ~ 3, without discussing possible 
evolution between the two redshifts. In this paper, we show 
the results of our stellar population analysis of these LAEs. 
We also compare the stellar populations of LAEs with those 
of other high redshift galaxy populations. The contribution 
of LAEs to the cosmic star formation density and the stellar 
mass density is also evaluated. 

The outline of this paper is as follows. In Section 2, we 
present our data and LAE sample. The SED fitting method 
is described in Section 3. In Section 4, we present and discuss 
our SED fitting results. A summary is given in Section 5. 
Throughout this paper, we use magnitudes in the AB system 
and assume a flat universe with (Q m , ^a, h) = (0.3, 0.7, 0.7). 



2 DATA 

2.1 Optical and NIR Images 

Deep BVRi'z' images of the SXDF were taken 
with Suprime-Cam on the Subaru Telescope by 
the Subaru/XMM-iV e mfon Deep Sur v ey project 
jSXDS; Furusawa et all |200S( ). lOuchi et all (^OOSf ) com- 
bined this public data set with their own imaging data 
taken with Suprime-Cam through two narrowband filters, 
7VB503 (A c = 5029A, AA = 74A) and AB570 (A c = 5703A, 
AA = 69A), and constructed samples of 368 2 = 3.1 and 116 
2 = 3.7 LAEs over a sky area of ~ 1 deg 2 . These numbers 
are the logical sum of the spectr oscopically-c o nfirm ed and 
photometrically-selected LAEs in lOuchi et all |200Sf ). They 
are selected as objects with narrow-band excesses and 
IGM absorption features shortward of Lyman a seen in 
color-color diagrams. The line-of-sight depth corresponding 
to the FWHM of the narrow-band filter is 55.7 comoving 
Mpc for NB503 and 43.2 comoving Mpc for NB570. For all 
LAEs, 2"-diam etre aperture mag nitudes were measured in 
each bandpass (|Ouchi et al.1120081 ). 

About 77% of the SXDF Suprime-Cam field was im- 
aged in the J, H, and K bands with the wide-field near in- 
frared c amera WFC AM on th e UKIRT in the UKIDSS/UDS 
project (|Lawrence et alj|2007h . The UKIDSS/UDS is under- 
way, and we use Data Release 3 for this study. We align the 
J, H, K images with the SXDS optical images using com- 
mon, bright stars, and then smooth them with Gaussian 
filters so that the PSF sizes of the J, H, K images match 
those of the optical images (FWHM Ril."0). The 3<r limit- 
ing magnitudes over a 2"-diametre aperture are calculated 
to be mj = 24.5, m» = 24.2, and rriK = 24.1. Because the 



1 http : //ssc . spitzer . caltech . edu/legacy/abs/dunlop .html 
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Table 1. Summary of Stellar Population Analysis of LAEs 



Reference 


Field 


Area 
[aremm J 


rcdshift 


JV ' 


™t 

three, 
[ma 


hold 

s] 


Bands 


Remark 


This Study 


SXDF 


2340 


3.1 


5(200) 


K (3ct) = 


= 24.1 


K, i, z, J, H, K, [3.6], [4.5], [5.8], [8.0] 


1 


Gawiser et al. (2006) 


ECDF-S 


992 


3.1 


(18) 






U, B, V, R, I, z, J, K 


2 


Gawiser et al. (2007) 


ECDF-S 


992 


3.1 


(52) 


[3.6](2<x) 


= 25.2 


U,B,V,R,I,z,J,K,[3.6], [4.5], [5.8], [8.0] 


3 


Lai et al. (2008) 


ECDF-S 


992 


3.1 


18(76) 


[3.6](2<r) 


= 25.2 


B, V, R, I, z,[3.6], [4.5] 


4 


Nilsson et al. (2007) 


GOODS-S 


44.4 


3.15 


(23) 


K„(3<t) 


= 23.4 


U, B,V,i,z', J, H,K S , [3.6], [4.5], [5.8], [8.0] 


5 


This Study 


SXDF 


2340 


3.7 


6(61) 


K (3<t) = 


= 24.1 


R, i, z, J, H, Jf,[3.6], [4.5], [5.8], [8.0] 


1 


Pentericci et al. (2009) 


GOODS-S 


160 


3.5 - 6 


((38)) 








(i 


Finkelstein et al. (2009) 


GOODS-S 


Kill 


4.5 


11(3) 






NB, V, i' ,z', J, H, if ,[3.6], [4.5] 


7 


Pirzkal et al. (20071 


HUDF 


11 


4.00 - 5.76 


5(4) 






z, J, H, K s ,[3.6], [4.5] 


8 


Lai et al. (2007) 


GOODS-N 


160 


5.7 


3 


[3.6](3a) 


= 26.2 


i', z',[3.6], [4.5] 


9 


Chary et al. (2005) 


Abcll 370 




6.56 


1 








10 



' Number of objects detected in at least one rest-frame optical band. The number in the parentheses is the number of objects not 
detected in any of the rest-frame optical bands, 
t Rest-frame optical detection threshold for LAEs, if given. 
* Bands used for SED fitting, if given. 
A"-detected objects are treated individually, while i<"-undetected objects are stacked. 

2 All objects are stacked. 

3 All objects are stacked. 

4 IRAC-detected and IRAC-undetected objects are stacked separately. 

5 All objects are stacked. 

6 LBGs with EW(Lya) > 20 A, and all are treated individually irrespective of the NIR. magnitude. 

7 All objects are treated individually irrespective of the NIR magnitude. 

8 All objects are treated individually irrespective of the NIR magnitude. 

9 All objects are treated individually. 

10 HCM 6A, a lensed LAE, detected at 3.6 and 4.5/jm. SED fitting is performed considering the effect of Ho+ [Nil] emission in 
the 4.5/^m band. 



zero-point magnitudes for the J, H, K images are given in 
the Vega system, we convert them int o AB magnitudes us - 
ing the offset values given in Table 7 of lHewett et al.l (|2006h . 

The SpUDS covers 0.65 deg 2 of the overlapping area of 
the SXDS and UDS fields (Figure [J). This 0.65 deg 2 area 
corresponds to an effective survey volume of 4.6 x 10 5 Mpc 3 
for z = 3.1 LAEs and 4.1 x 10 5 Mpc 3 for z = 3.7 LAEs, re- 
spectively. All of the SpUDS IRAC images are geometrically 
matched to the optical images. We calculate the 3u limiting 
magnitude over a 3"-diametre aperture to be rri3.6 = 24.8, 
7T14.5 = 24.5, m 5 .8 = 22.7, and ms.o = 22.6 in the 3.6, 4.5, 
5.8, and 8.0/im IRAC bands, respectively. 



2.2 Photometry 

In this paper, we only analyse LAEs in the overlapping area 
of 0.65 deg 2 where the Suprime-Cam, WFCAM, and IRAC 
data (either channel 1+3 or 2+4) are all available. We per- 
form J, H, and K photometry with a 2"-diametre aperture 
at the position of LAEs in the narrowband images, using the 
IRAF task apphot. 

We then convert 2"-diametre aperture magnitudes in 
the optical and JHK bands into total magnitudes in the 
following manner. First, we select 20 bright and isolated 
point sources in the i?-band image, the deepest among all 
the images, and measure fluxes in a 2" aperture and in a se- 
ries of larger apertures up to 5" with an interval of of 0.1". 
Since we find the fluxes to level off for > 4" apertures, we 
define 4"-aperture magnitudes as total magnitudes. Then, 
we select 100 point sources in the 7i-band image, measure 
fluxes over 2" and 4" apertures, and calculate an accurate 
offset between these two aperture magnitudes to be 0.17 
mag. For each of the optical and JHK bandpasses, we sub- 
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Figure 1. Observed field in the WCS coordinate system. The 
blue cross-shaped region is the SXDF imaged with Suprime-Cam. 
The green square corresponds to the UKIDSS/UDS (WFCAM) 
field. The red solid line and red dashed line outline the SpUDS 
fields with IRAC 3.6/^m and 5.8/^m imaging and 4.5/^m and 8.0^tm 
imaging, respectively. 



tract 0.17 mag from 2"-aperture magnitudes to obtain total 
magnitudes for our LAEs. 

For the Spitzer/IRAC four bands, we measure 3"- 
diametre aperture magnitudes for each LAE and converted 
them to total magnitudes by applying the aperture correc- 
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tion given by MUSYC survej0. The correction values are 
0.52, 0.55, 0.74, and 0.86 mag for 3.6/xm, 4.5/zm, 5.8zjm, 
and 8.0/xm, respectively. 

2.3 The LAE samples used for the stellar 
population analysis 

Table [5] summarises our LAE samples used for the stel- 
lar population analysis. Our original sample consists of 224 
LAEs at z = 3.1 and 78 LAEs at z = 3.7, among which 
35 and 21 have been spectroscopically confirmed, including 
two newly confirmed objects using Magellan/IMACS (Sec- 
tion [23). 

We regard an object as LAE with AGN and exclude it 
from our analysis if it has a counterpart eith er in the XMM- 
Newton X-ray c atalogue dUeda et al.| [2008') or in the VLA 
radio catalogue l|Simpson et al.ll2006h . or its spectrum has 
emission lines typical of AGN, such as like SilV, CIV, and 
Hell (for details, see Section 4.6 in lOuchi et~afl 120081 ') . It 
should be noted, however, that we cannot completely ex- 
clude contamination from X-ray-faint, radio-quiet, and/or 
heavily obscured AGNs. 

In the stellar population analysis, we treat objects 
detected in the K band individually, while we stack K- 
undetected objects at each redshift to make an average SED 
and fit it with stellar population synthesis models. The K 
band corresponds to the rest-frame V and B bands for 
2 = 3.1 and z = 3.7, respectively, both of which are red- 
ward of the 4000 A break. We find that the stacked objects 
of the two redshifts are both detected in K. Since as high 
as 96% of our sample are undetected in K (see below), the 
stacked SEDs of if-undetected objects is considered to rep- 
resent closely the whole LAE population at the two redshifts. 

2.4 iC-detected LAEs 

Among a total of 224 (78) LAEs at z = 3.1 (3.7), 9 (9) are 
found to be brighter than the if-band 3er magnitude (i.e., 
24.1 mag). Among these 18, 2 (1) objects at z = 3.1 (3.7) are 
found to host AGN (judged from the X-ray and radio data 
and spectroscopic data as described above), and 2 objects 
at 2 = 3.7 are found to be significantly confused by their 
neighbouring objects, we do not use them for the stellar 
population analysis. We then visually inspect the images of 
the remaining 13 if-bright objects, and find that one object 
is likely to be a false detection, because it is not detected in 
any of the other IR bands (J, H, and IRAC bands) and it 
falls on a low quality part of the if-band image. Next, we 
examine the SEDs of the remaining 12 objects over the range 
from the B band to the IRAC 8.0^im band, to find that one 
has an SED consistent with a late-type star, with a peak 
near the K band. After removal of these two objects, the 
number of LAEs brighter than K(3a) is five for z = 3.1 and 
six for 2 = 3.7. Hereafter, we call them "if -detected LAEs". 
Table [3] summarises their photometry, L(Lya), EW(Lya), 
and redshift. 

We have found that the source center of NB503-N-21105 
in the narrow band (Lya) offsets from those in the broad 

2 http : //www . astro . yale . edu/dokkum/SIMPLE/ 
release.! . l/doc/00README_photometry 



bands (continuum) by about one arcsec. We thus cannot 
rule out the possibility that NB503-N-21105 is an LAE with 
a chance projection of a galaxy along the line-of-sight. How- 
ever, the two components are more likely to be physically 
connected, because the narrow-band image is elongated to 
the direction of the offset. 

NB570-N-32295 has a power-law like spectrum over the 
whole wavelength range up to 8.0/^m. Although we find that 
this spectrum is reproduced reasonably well by a stellar sys- 
tem (see Section T3.2p . we cannot completely rule out a pos- 
sible contamination from an AGN. 



2.5 New Spectroscopy of Two Red Objects 

Six out of the 11 i f -detected LAEs ha ve been spectroscop- 
ically identified by lOuchi et al.l (|200&t ) and Ouchi et al. in 
preparation. Among the remaining five, we selected two red 
objects, NB503-N-21105 and NB570-N-32295, and carried 
out spectroscopy with the Inamori Magellan Areal Camera 
and Spectrograph (IMACS) on the Magellan I Baade 6.5m 
telescope at the Las Campanas Observatory. The observa- 
tion was made on 2008 December 19 with the WB4800 - 
7800 filter and the 300 lines mm -1 grism whose blaze angle 
is 17.5 degrees. We chose a slit width of 1.0". The on-source 
exposure time was 15300 seconds under the 0.5" - 0.75" 
seeing condition. The spectral coverage was 4800 - 7800 A. 
The spectral resolution and the corresponding velocity res- 
olution are R ~ 890 and Av ~ 340 km s" 1 at 5029 A, and 
~ 1000 and 300 km s -1 at 5703 A. The data were reduced 
with the COSMOS pipeline version 2.12^1 

We present the reduced 2D and ID spectra in Figure [2] 
Both objects have a single, strong emission line (> 6a de- 
tection for NB503-N-21105 and > 4a detection for NB570- 
N-32295) and have no other emission line or detectable con- 
tinuum emissiorQ. The central wavelength of the single line 
is 5035.4 A for NB503-N-21105, and 5693.9 A for NB570- 
N-32295. We conclude that both objects are a real LAE 
from the following discussion. If NB503-N-21105 is an [O III] 
([Oil]) emitter, Ha ([Olll]) emission should be detected 
at ~ 6629 A (6759 A). Similarly, if NB570-N-32295 is an 
[Olll] ([Oil]) emitter, Ha ([Olll]) emission should be seen 
at ~ 7482 A (7661 A). We summarise the spectroscopic re- 
sults in Table g] 



2.6 Stacking of if-undetected LAEs 

Among the LAEs fainter than the K(3a) magnitude, 200 
(61) at 2 = 3.1 (3.7) have all 10 broadband data from 
R to 8.0/im and do not show AGN features, to be used 
for the stellar population analysis. We call these LAEs 
"if- undetected LAEs". We make median-stacked multi- 
waveband images separately for the two redshifts. Table [5] 
summarises their photometry, average L(Lya), EW(Lya) 



li http : //www. ociw. edu/Code/cosmos 

4 There is an emission-like feature around [C IV] in the spectrum 
of NB570-N-32295. We find, however, that there is a strong sky 
emission line at this wavelength and that this feature is statis- 
tically insignificant, being consistent with a residual of sky sub- 
traction. 
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Table 2. Summary of our samples used for the stellar population analysis 



LAE sample 


A"- band 


N 


with 


referred to as 




detection 




spec-2 




2 = 3.1 


Total 




224 


35 




LAE without AGN 


Yes 


5 


5 


If -detected LAEs 


LAE with AGN 


Yes 


2 


2 




noise (false detection) 


Yes 


1 







late-type star 


Yes 


1 







LAE in the SpUDS fileds and without AGN 


No 


200 


27 


^-undetected LAEs 


LAE out of the SpUDS fileds and/or with AGN 


No 


15 


1 




2 = 3.7 


Total 




78 


21 




LAE without AGN 


Yes 


6 


3 


if -detected LAEs 


LAE with AGN 


Yes 


1 


1 




confusion 


Yes 


2 


2 




LAE in the SpUDS fileds and without AGN 


No 


61 


11 


^-undetected LAEs 


LAE out of the SpUDS fileds and/or with AGN 


No 


8 


1 





Table 3. If-detected LAE sample 



Object Name 


R 


i 


z 


.7 


II 


K 


3.6^tm 


4.5/tm 


5.8^tm 


8.0(im 


log L(hya) 
[erg s _1 ] 


EW(Lya) 

[A] 


2 


Ref. 


NB503-N-21105 (Rl) 


26.2 


25.9 


26.0 


24.3 


23.6 


23.9 


22.4 


22.5 


22.8+ 


22.9+ 


42.49 


73 


3.142 


(1) 


NB503-N-42377 


23.8 


23.8 


23.9 


23.8 


23.8+ 


23.0 


23.7 


24.2+ 


23.1 + 


22.3+ 


43.50 


77 


3.154 


(2) 


NB503-S- 45244 


24.6 


24.6 


24.8 


99.9+ 


25.1 + 


23.6 


24.7+ 




99.9+ 




43.21 


86 


3.156 


(2) 


NB503-S-65716 


23.9 


23.9 


24.1 


24.3 


23.8 


23.5 


24.1 




99.9+ 




43.38 


65 


3.114 


(3) 


NB503-S-94275 


24.3 


24.3 


24.5 


24.5+ 


24.8+ 


23.6 


24.9+ 


24.4+ 


99.9+ 


23.7+ 


43.56 


147 


3.102 


(2) 


NB570-N-32295 (R2) 


27.1 


26.5 


26.2 


25.3+ 


24.6+ 


23.5 


22.6 


22.3 


21.6 


21.5 


42.74 


225 


3.684 


(1) 


NB570-S-84321 


24.9 


24.8 


25.0 


24.4+ 


24.0 


23.7 


24.3+ 


23.7 


26.7+ 


24.2 + 


42.99 


55 


3.648 


(2) 


NB570-W-55371 


24.6 


24.6 


24.6 


25.3+ 


24.1 + 


23.3 


23.8 


23.9 


23.2+ 


99.9+ 


42.71 


20 


3.699 


(2) 


NB570-C- 24119 


25.5 


25.2 


25.1 


26.8+ 


24.7+ 


23.2 


23.0 


22.5 


22.0+ 


22.0+ 


42.58 


34 


3.69* 




NB570-S-88963 (R3) 


26.1 


25.4 


25.4 


24.2 


25.2+ 


23.2 


23.1 


22.6 


22.8+ 


22.8+ 


42.83 


106 


3.69* 




NB570-W-59558 (R4) 


25.7 


25.4 


25.2 


24.2 


24.0 


23.0 


23.3 


23.1 


99.9+ 


23.8+ 


43.12 


150 


3.69* 





NOTES: All magnitudes are total magnitudes. No value means that the object is out of the image. 
Negative fluxes ha ve been replaced with 99.9 mag. Sources of spectroscopic redshifts are (1) this study, (2) 
lOuchi et al. | l l200gt) , (3) Ouchi ct al. in preparation. 
+ Fainter than the 3c limiting magnitude. 

*Not spectroscopically confirmed, we give z = 3.69, which corresponds to the central wavelength of NB570. 



Table 4. Spectroscopic results of two red LAEs 



Object Name 


RA (J2000) 


Dec (J2000) 


FWHM(Lya)t 


z 




[h:m:s] 


[d:m:s] 


[km s" 1 ] 




NB503-N-21105 (Rl) 


2:18:42.186 


-4:46:38.54 


205 ± 129 


3.142 


NB570-N-32295 (R2) 


2:17:25.630 


-4:44:33.57 


629 ± 201 


3.684 



+Aftcr correction for the instrumental broadening on the assumption of a Gaussian profile. 



Table 5. Stacked median LAE sample. 



Object Name 


R 


i 


z 


J 


H 


K 


3.6^tm 


4.5^tm 


5.8/^m 


8.0^m 


log £(Lya) 


EW(Lya) 


























[erg s _1 ] 


[A] 




NB503-if-undetected 


26.9 


26.9 


27.1 


27.4+ 


28.8+ 


26.5 


27.1 


27.4+ 


99.9+ 


99.9+ 


42.54 


155 


3.14* 


NB570--fS"-undetected 


26.4 


26.3 


26.4 


26.9+ 


27.2+ 


25.4 


26.5 


29.2+ 


25.5 + 


99.9+ 


42.82 


135 


3.69* 



NOTES: All magnitudes are total magnitudes. Negative fluxes have been replaced with 99.9 mag. 
+Fainter than the 3<r limiting magnitude. 

* Redshifts corresponding to the central wavelengths of NB503 and NB570. 
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Figure 2. Spectra of NB503-N-21105 (top) and NB570-N-32295 (bottom). For each object, the top and middle panels show the two- 
and one-dimensional spectra, respectively. The one-dimensional spectra have been smoothed with a 3 pixel boxcar. The dashed lines 
with the legend correspond to the wavelengths of typical emission lines from AGN. Inserted is a zoom up around Lyo line. The bottom 
panel shows the sky background with an arbitrary normalisation. 
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lar p opulations by the standard SEP fitting method 
(e.g., iFurusawa et all l200d ; iPapovich et al.l 120011 ). We use 
the stellar population synthesis model of GALAXEV 
|BruzuaI fc Cha riot 2003, hereafter BC03) to produce model 
SEDs. Most of the previous studies have used BC03 (Ta- 
ble [SJ. We make a large set of mass-normalised model 
SEDs, varying star-formation timescale, age, and dust ex- 
tinction (or E(B — V)*). These SEDs are then redshifted 
to z = 3.1 and z = 3.7 and convolved with ten bandpasses, 
R, i, z, J, H, K, and IRAC four bands, to calculate flux den- 
sities. For each LAE, we search for the best-fit SED that 
minimizes 



Figure 3. Normalised histograms of the rest-frame Lya EWs for 
all (black) and A"-undctectcd (grey) LAEs at z = 3.1 (left) and 
3.7 (right). 



[/gi - AWgLi (age. E{B - V%r)] ' 

* 2 (/iS) 



(1) 



derived from narrow-band and i?-band magnitudes, and red- 
shift. In Figure [3] we show the histograms of the rest-frame 
Lya EWs for all and /("-undetected LAEs, to confirm that 
the ^-undetected LAEs are representative of the whole sam- 
ple. 



2.7 Applications of Broadband Color Selections 

2.7.1 LBG Selection 

Wc examine whether the six _R"-detected LAEs at z = 3.7 are 
selected as LBGs, by appl ying the BRi ' colour selection cri- 
teria for z ~ 4 LBGs (e.g. lOuchi et ai]|2004l : lYoshida et~aH 
I2006T I. None of them are found to meet the criteria. This is 
probably because our objects tend to have redder UV slopes 
than typical LBGs, or the B-band image is not deep enough 
to detect Lyman break, or the redshift of z — 3.7 is close to 
the edge of the redshift range covered by the BRi' selection 
criteria. 

In addition, we apply the BRi' colour criteria to the 
61 ^-undetected LAEs at z = 3.7 and find that only seven 
satisfy the criteria. We also find that the stacked LAE at 
2 = 3.7 does not meet the criteria. 

2.7.2 DRG Selection 

We then apply the colour criterion for distant-red galax- 
ies (DRGs), J — K > 1.3, to our sample, and find that 
none of the if-detected LAEs and the two stacked LAEs 
satisfies the DRG criterion. We also apply the colour cri- 
terion to the best-fit model SEDs for these LAEs derived 
from our SED fitting (Section[3]). Among the 11 7f -detected 
LAEs, model SEDs for three (NB503-N-21105, NB570-N- 
32295, and NB570-C-24119) meet the criterion. The former 
two have very red SEDs and have very large dust extinction 
(E(B - V)* ~ 0.65). We will focus on this kind of LAEs in 
Section [O] 



3 SED FITTING 
3.1 Method 

After obtaining rest-frame UV-to-optical photometry for 
the if-detected and stacked LAEs, we analyse their stel- 



where 



is the observed flux density in the i-th bandpass, 



/model i s the mass-normalised model flux density in the i-th 
bandpass, r is the star- formation timescale, and <r(/ffl) is 
the photometric error in the i-th bandpass. Because M a t ar is 
the amplitude of a model SED, we obtain the best-fit M atar 
for each set of (age, E(B — V) + , r) by solving dx 2 /9M BtSlT = 
0, and calculate \ 2 ■ Then we search for the set of the best- 
fit parameters that gives the minimum \ 2 ■ The errors in 
the best-fit SED parameters correspond to la confidence 
interval (Ax 2 < 1) for each parameter. 

We do not use either B- or V-band photometry in 
the fitting, since _B-band photometry suffers from the IGM 
absorption shortward of the Lyman a wavelength and V- 
band photometry is contaminated from Lya emission. The 
amount of the IGM absorption considerably differs by the 
line of sight. Lya fluxes can be estimated from narrow and 
broad band photometry, but they have an unc ertainty of a 
factor of ~ 2 (see Figure 15 of lOuchi et al]l2008h . The short- 
est broadband we use for SED fitting is the R band, which 
corresponds to 1585 A and 1380 A for z = 3.1 and 3.7 LAEs. 
These wavelengths are very close to Lya and short enough 
to cover the wide range of SEDs. 

W e adopt Salpeter's initial mass function (jSalpeterl 
Il955h with lower and upper mass cutoffs of 0.1 and IQQMq. 
We fix the metallicity to Z/Zq — 0.2, considering the fact 
that LBGs at z ~ 3 tend to have subsolar metallicies 
at Z/Z Q ~ (1/4) - (1/3) (a strongly lensed LBG, cB58: 
iPettini et aHl200Cl : iTeplitz et alJl200Ch and Z/Z G) ~ (1/10) 
- (2/3) (four bright LBGs: IPettmi et alj|200lh . Although 
IShaplev et al.l (|2004h reported that very massive z ~ 2 LBGs 
with M B tar ~ 10 11 Ma have approximately the solar metal- 
licity, UrbieFal] (|2006l ) found that less massive ( < 1O 1O M ) 
LBG s at similar redshift s have subsolar metallicities (see 
also. iMaiolino et al.ll2008T ). It appears to be reasonable that 
the metallicites of typical LAEs are lower than (and at most 
comparable to) those of typical LBGs. 

We also examine models with Z/Zq — 0.005 and 1, 
and find that while the best-fit stellar mass is insensitive 
to metallicity over this range, the best-fit age and dust ex- 
tinction become higher with decreasing metallicity. However, 
the dependencies of age and dust extinction on metallicity 
are not so strong: for both parameters, the change from the 
best-fit value for the Z/Zq = 0.02 model is within its la 
errors. 
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Table 6. Summary of Stellar Population Analysis of LAEs: Assumption on SED Modeling 



Reference 


Modcl^' 


IMF 


SFH* 


Metallicity 


Extinction* 


Remarks 












Curve 




This Study 


BC03, CB08 


Salpctcr 


const, exp 


0.005, 0.02, 0.2, 1 


COO 


(O09) 


Gawiser et al. (20061 


BC03 


Salpctcr 


const 


1 


C97 




Gawiser et al. (2007) 


M05 


Salpctcr 


two-burst 


0.02 - 1 


COO 




Lai et al. (2008) 


BC03 


Salpctcr 


const 


1 


COO 




Nilsson et al. (2007) 


BC03 


Salpctcr 


const 


0.005, 0.2, 1 


CF00 




Pentericci et al. (2009) 


BC03, M05, CB08 


Salpctcr 


exp 


0.02, 0.2, 1, 2.5 


COO 


(P09) 


Finkelstein et al. (2009) 


BC03 


Salpctcr 


cxp, two bursts 


0.005, 0.02, 0.2, 0.4, 1 


C'!)4 


(F09) 


Pirzkal et al. (2007) 


BC03 


Salpctcr 


ssp, cxp, 2bp 


0.005 - 2.5 


COO 




Lai et al. (2007) 


BC03 


Salpctcr 


ssp, const 


0.005, 1 


coo 




Charv et al. (2005) 


BC03 


Salpctcr 




0.02 







t Population synthesis model. "BC03", "M05", and "CB08" represent iBruzual fc Chariot! [120031) , iMarastonl j2005h . and the 
unpublished Chariot & Bruzual model, respectively. 

^ Star formation history, "ssp": instantaneous starburst, "const": constant star formation history, "exp": exponentially 
decaying star formation history (and r is e-folding time), "two-burst": superposition of an old instantaneous burst component 
and a young "exp" component, "two bursts": superposition of an old "exp" component (age = 1.4 Gyr and r = 10 s yr) and a 
young "exp" component (r = 10 5 yr), "2bp": superposition of two different instantaneous starbursts. The e-folding times 
examined are: r = 1, 10, 10 2 , 10 3 Myr (O09), r = 0.1, 0.3, 0.6, 1, 1.5, 2, 3, 4, 5, 7, 9, 13, 15 Gyr (P09), r = 10 5 , 10 6 , 10 7 , 
10 8 , 4 X 10 9 yr (F09). 

* Dus t extinction law. COO: ICalzetti et alj <200CD . C97: ICalzetti et alj jl997l) , C94: ICalzetti et all [119941) , CF00: ICharlotfcFalll 

feoodi 



The star formation timescales examined are (i) constant 
star formation and (ii) exponentially decaying star forma- 
tion with four e-folding times of r = 0.001, 0.01, 0.1, 1 Gyr, 
and we search for the best-fit SED separately for the cases (i) 
and (ii). As explained in the next section, we adopt the re- 
sults for constant star formation for our discussion of stellar 
populatio ns. For dust extinct ion, we use Calzetti's extinc- 
tion law (|Calzetti et al.ll2000l ) and vary E(B - V)* over 
and 1.50 with an interval of 0.01. 

Models of constant star formation have three free pa- 
rameters, stellar mass, age, and dust extinction, while mod- 
els of exponentially decaying star formation have one more 
parameter, e-folding time. 



3.2 Results 

Figures [4] and [5] show the results of the SED fitting for ob- 
jects at z = 3.1 and 3.7. For each object, the blue and red 
curves correspond respectively to the best-fit SEDs for con- 
stant star formation and for exponentially decaying star for- 
mation. Both SEDs give similarly good fits for most of the 
objects, implying a difficulty in constraining the star forma- 
tion history. 

The best-fit parameters are summarised in Table [7] for 
constant star formation and Table[8]for exponentially decay- 
ing star formation. The reduced \ 2 tends to be slightly larger 
for the model of exponentially decaying star formation. This 
may be because the e-folding time interval adopted is too 
large (1 dex) to fine-tune model SEDs. In any case, adding 
the e-folding time as a free parameter does not significantly 
improve the fit. 

It is found from these tables that the two models give 
very similar stellar masses, and thus that the stellar mass 
is a robust parameter against the assumed star formation 
history. Dust extinction is also relatively insensitive to the 
star formation history. The only exception is NB570-W- 
59558, for which E(B — V)* = 0.63 for constant star for- 
mation and 0.00 for exponentially decaying star formation. 



On the other hand, age and SFR largely differ between 
the two models for some objects. The ages of two objects 
(NB570-C-24119 and NB570-W-59558) differ by more than 
one order of magnitude, and the SFRs of seven objects 
(NB503-S-65716, NB503-S-94275,NB503-N-21105, NB503- 
A'-undetected, NB570-C-24119, NB570-N-32295, NB570-W- 
59558) differ by more than one order of magnitude. All but 
one of these seven objects have r = 0.001 Gyr, i.e., burst- 
like star formation histories. In other words, these objects 
are equally well fitted by two extremes, constant star for- 
mation and burst-like star formation, and our data cannot 
discriminate betw een them. It will be resolv ed by deeper 
NIR imaging (e.g.. |Pozzetti fc Mannucclll2000l ) and/or other 
independent observations such as submillimetre imaging. 

In the following discussion, we will adopt the results for 
constant star formation, because (i) the values of reduced 
X 2 are lower for constant star formation for all but one ob- 
ject and (ii) it is easy to compare our results with previous 
SED studies on LAEs since most of them have also adopted 
constant star formation. 

Most of the LAEs have a brighter A-band magnitude 
than the best-fit SED. The reason for this discrepancy is not 
clear, but a possible reason is that the A-band flux density 
is affected by [O III] and H/3 emission. On the basis of the 
relations between SFR and [O Hi] and LL 3 emission strengths 
|KennicuttJll998l : iMoustakas et alj|2~006f ) , we find that these 
lines may significantly increase A-band flux densities if SFR 
is higher than ~ 10 3 [Mg yr -1 ]. However, since most of the 
objects have SFR < 1O 3 M0 yr -1 , this effect cannot explain 
the observed excesses. Another possible reason is that a rel- 
atively faint magnitude, K(3a) = 24.1, is adopted for the 
boundary of i\-detection. Our A-detected sample may in- 
clude objects whose K magnitude happens to be brighter 
than the true value due to positive sky noise at the posi- 
tion of the objects. However, most of the i\-detected ob- 
jects are also detected in at least one of the J, H, and the 
four IRAC bands, suggesting that they are truly bright in 
the near-infrared wavelengths and thus that the systematic 
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Figure 4. Best-fit model SEDs (curves) and observed photometry (squares) for X-detected LAEs. The blue and red curves correspond 
to constant star formation and exponentially decaying star formation, respectively. Data shown by open squares are not used for the 
SED fitting. 



brightening of K is not very strong. We infer that the sys- 
tematic brightening of K, if any, does not significantly af- 
fect the best-fit SEDs which are determined from the overall 
shapes of observed SEDs from R to IRAC 8.0/im bands. In- 
deed, we perform the SED fitting excluding K photometry, 
and find that for any object the best-fit parameters which 
are different from the original values only within the la er- 
rors, although the best-fit stellar masses tend to be slightly 
(~ 0.05 dex) smaller. 

Most of the errors in the best-fit parameters of NB503- 
if-undetected are larger than those of NB570-A"-undetected 
in spite of a larger number of objects being stacked in 
NB503. This is because NB503-7f -undetected is fainter in 
most bandpasses as shown in Table 5. The z = 3.1 LAE sam- 
ple includes fainter objects than the z — 3.7 sample since the 
limiting m agnitude in NB503 is fainter by 0.6 mag (see Sec- 
tion 2.2 in lOuchi etalll2008h . Thus, the /("-undetected ob- 



jects in the z = 3.1 sample are on average fainter than those 
in the z = 3.7 sample. We also note that the reduced \ 2 is 
smaller for NB503-A"-undetected. We infer that this is not 
only because the observed A"-band flux density for NB570- 
A-undetected is much brighter than the best-fit SED (see 
Figure [5]), but also because the magnitude errors of NB503- 
A"-undetected are larger. 

Some new population synthesis models include ther- 
mally pulsating asymp t otic giant bran ch (TP-AGB) stars 
(e.g., iMarastonl 12003 ; iBruzuall I2007T I. Inclusion of TP- 
AGB stars should have little effect on the SED fitting 
for most LAEs, because LAEs are in general very young 
and have sub solar metallicities (e.g.. IPentericci et al J 120091 : 
lGawiserll2009l ). In our sample, however, four objects (NB503- 
N-42377, NB503-S-45244, NB570-W-55371, and NB570-C- 
24119) have relatively old ages of > 300 Myr and thus their 
NIR fluxes could be dominated by TP-AGB stars. We fit 
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Figure 5. Same as Figure [4] but for iC-undetected (stacked) LAEs. 



Table 7. SED Fitting Result (Z = 0.2Z Q , constant SFH) 



Object Name 


log M s tar 


E{B — V)* 


log(Age) 


log(SFR) 


xl* 






[mag] 


[yr] 


[M @ yr- 1 ] 




z = 3.1 


NB503-N-21105 (Rl) 


io.43lg;S 


u '' u -0.05 




o yc-4-0.34 
— .40 


2.36 


NB503-N-42377 


9.82lH 9 o 


n nn+ 002 
u - uu -o.oo 


8.51 j 


0.10 
0.20 




2.75 


NB503-S-45244 


9.58±g : g 


n nn+ 007 
u - uu -o.oo 


o fil+0.25 
8 - bi -0.70 


i nfi +0 - 36 

1-UD_q Q1 


2.17 


NB503-S-65716 




n in+0.05 
u - iJ -0.03 


6.98j 


0.40 
0.34 


9 99 +0.46 
z - zz -0.26 


1.91 


NB503-S-94275 


8.97±°il 


n n + 07 
u - ii -o.n 


7.10j 


1.26 
0.48 


1 HS+ - 59 


0.96 


NB503-.ff-undetected 


8 19+0.26 


n 03+ ' 05 


7 O1+0.45 
'•81-0.81 


0.36i°;f 7 


1.39 



3.7 



NB570-N-32295 (R2) 
NB570-S-84321 
NB570-W-55371 
NB570-C-24119 
NB570-S-88963 (R3) 
NB570-W-59558 (R4) 
NB570-if-undetected 



9.59lg;f 6 
10.19±g;i? 

7+0.17 
'-0.15 
1U.O^_ 33 
,+0.03 
-0.69 



10.67"* 



10.71] 



.50 



+0.20 
-0.10 



67+0- 04 
u - D '-0.25 
HO. 05 
-0.08 
4+0.04 
-0.04 
HO. 07 
-0.08 
}+0.06 
-0.10 

3+0.01 

-0.16 

^+0.04 
-0.03 



0.32Z 
0.04j 
0.29^ 
0.62^ 
0.63^ 
0.19j 



7.40 
6.68 
8.96 
8.71 
6.56 
6.16 
6.76 



+ 1.81 
0.38 
o.sr, 
0.12 
+0.25 
0.35 
+0.50 
0.50 
+0.44 
0.08 
+0.74 
1.06 
+0.16 
0.22 



3.23 
2.91 
1.33 
2.06 
3.98 
4.55 
1.74 



+0.27 
1.19 
+0.35 
0.73 
+0.19 
0.15 
+0.33 
0.33 
+0.37 
0.75 
1.09 
1.42 
+0.42 
0.22 



0.57 

1.34 
1.15 
1.39 
1.93 
2.14 
3.27 



'Reduced \ squares. The degree of freedom is 5 (= 
7 (= 10 - 3) for the others. 



8-3) for NB503-S-45244 and NB503-S-65716, and 



a new version of GALAXEV which includes TP-AGB stars 
(Chariot & Bruzual, in preparation) to these four, and ob- 
tain almost the same results as for BC03, except decreases 
in stellar mass of 0.2 - 0.3 dex. 



4 RESULTS AND DISCUSSION 

Figure [6] shows the distributions of the best-fit parameters 
derived from the SED fitting. The stellar masses of the K- 
detected LAEs are distributed over M star ~ 10 9 - 10 10,5 M Q , 
while those of the JL-undetected LAEs are much less massive 



10 8 - 1O 8 ' 5 M ). The ^-detected LAEs have a large 



variety of dust extinction with < E(B — V)* < 0.7, while 
the /^-undetected LAEs have mild dust extinction. The ages 
of our LAEs are distributed over around 10 6 - 10 9 yr. The 
age distribution of the /('-detected objects could be bimodal, 
although its statistical significance is low. The ages of the K- 
undetected LAEs are around the median of the /^-detected 
LAEs. The SFRs of the Jf-detected LAEs span a wide range 
of 10 - 10 4 M o yr -1 , while those of the if-undetected LAEs 
are between 1 - lOOM© yr - . 

The E(B - V)* of NB570-if-undetected is significantly 
larger than that of NB503-if-undetected. This may merely 
reflect the fact that NB570-if-undetected is more massive 
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Table 8. SED Fitting Result (Z = O.2Z , exponentially decaying SFH) 



Object Name 


log M s tar 


E(B- V)* 


log(Age) 


log(SFR) 


e-folding time 


Xr 




[Mq] 


[mag] 


[yr] 


[M e yr" 1 ] 


[Gyr] 




2 = 3.1 


NB503-N-21105 (Rl) 


in is+ - 51 


64+ ' 11 






0.001 


2.71 


NB503-N-42377 


a -° z -0.12 


00+ 02 

u - uu -o.oo 




-, OO+0.09 

l-3»_o. 10 


1 


3.21 


NB503-S-45244 




n on+ 007 
u - uu -o.oo 


8.5lt°;f 7 




1 


2.72 


NB503-S-65716 




n n + 03 

u - ii -0.04 
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NB570-N-32295 (R2) 10.84to.37 0.48±o; 11 7.81±o 73 -22.9ll^' 29 0.001 0.56 

NB570-S-84321 ^-^to.to °- 31 to.30 6 - 68 ±o'.l6 2 - 75 -°J' 1 °- 01 L55 

NB570-W-55371 10.05±g;g| 0.00±{j;{$ 8.46±gg° 0.93±g;Jo 0.1 1.08 

NB570-C-24119 10.39±o.o7 °-29±o.06 7 -60to.il -12.92tg 5 4 17 0.001 1.53 
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NB570-W-59558 (R4) 10.46±g;°| 0.00t o ;°J 8.56i ;°g -12.08l^ 65 0.01 1.85 

NB570-i\-undetected 8.49lg;^ 0.19i 'og 6.74t°H 1.74±^|| 0.1 3.82 



'Reduced \ squares. The degree of freedom is 4(= 
6 (= 10 - 4) for the others. 
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Figure 6. Distributions of the best-fit parameters derived from the SED fitting. The shaded histograms are for the i\-detected LAEs 
at z = 3.1, and the open histograms are for the sum of the i\-detected LAEs at z = 3.1 and 3.7. Grey and black arrows indicate the 
/■("-undetected LAEs at z = 3.1 and 3.7, respectively. 



and thus has produced more metal. On the other hand, this 
could reflect some evolution of the LAE population between 
the two redshifts. Indeed, it is found in Table [5] that LAEs 
at z > 3.7 can have E(B — V)* comparable to or larger than 
that of NB570-Js:-undetected. 

We find that four objects have extremely large ex- 
tinction E(B - V)* ~ 0.65: NB503-N-21105, NB570-N- 
32295, NB570-S-88963, and NB570-W-5955© Their prop- 



erties (other than age) are clearly different from the other 
LAEs, as we can see in the following subsections. We label 
the four objects as Rl - R4 in most of the forthcoming fig- 
ures. These four red LAEs may be a distinct population, 
and we will discuss this in Section T4.2I 



5 This object is equally well fitted by two extremes: a young 
and dusty starburst galaxy with constant SFH, and an old and 
dead galaxy with exponentially decaying SFH. This is caused by 
the age-dust degeneracy. Although our data cannot discriminate 
between them, we treat this object as a young and dusty starburst 
galaxy in this section (see details in Scction l3.2t 



4.1 Correlations among Stellar Population 
Parameters 

Table|9]summarises the results of the stellar population anal- 
ysis of LAEs reported so far. In this subsection, we examine 
correlations among stellar population parameters, referring 
to the previous studies. 
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Table 9. Summary of Stellar Populations 



Reference 


z 




E(B - V)„ 


Age 


Z 


SFR 


T 


Remarks 






[10 9 Mg] 


[mag] 


[Myr] 


[Z Q ] 


[M yr _1 ] 


[Myr] 




This Study 


3.1 


0.13 


0.03 


65 


(0.2) 


2.3 


(oo) 


(O09-1) 


This Study 


3.1 


0.93 - 27 


0.00 - 0.70 


4.8 - 407 


(0.2) 


11 - 5.6 X 10 3 


(oo) 


(O09-2) 


Gawiser et al. (2(306) 


3.1 


0.5 


0.0 


90 


(1) 


6 


(oo) 




Gawiscr et al. ^007) 


3 . 1 


1 .0 


0.0 


20 


(1) 


2 


750 


(G07) 


Lai et al. (2008) 


3.1 


0.3 


0.0 


1.6 X 10 2 


(1) 


> 


(oo) 


(L08-1) 


Lai et al. (2008) 


3.1 


9 


0.0 


1.6 X 10 3 


(1) 


6 


(oo) 


(L08-2) 


Nilsson et al. (2007) 


3.1 


0.47 


0.07 


8.5 X 10 2 


(0.005) 


0.66 


(oo) 




This Study 


3.7 


0.32 


0.19 


5.8 


(0.2) 


55 


(oo) 


(O09-1) 


This Study 


3.7 


3.9 - 51 


0.04 - 0.67 


1.4 - 9.1 X 10 2 


(0.2) 


21 - 3.5 X 10 4 


(oo) 


(O09-2) 


Finkelstein et al. (2009) 


4.5 


0.084 - 6.1 


0.02 - 0.41 


2.5 - 500.0 


0.005 - 1 


1.30 - 5.78 


0.1 - 4 x 10 3 


(F09) 


Pirzkal et al. (2007) 


4.00 - 5.76 


0.007 - 1.8 


0.0 - 0.16 


1.0 - 20 


0.005 - 2.5 




1 - 12.5 


(P07) 


Lai et al. (2007) 


5.7 


17 - 39 


0.150 - 0.225 


720 - 900 


(0.005) 




(oo) 


(LOT) 


Chary et al. (2005) 


6.56 


0.84 


0.25 


5 


0.02 


140 




(C05) 



NOTES: Values in parentheses are the assumed values in the literature. 

REMARKS: (O09-1) i\-undetected LAEs. (O09-2) i\-detected LAEs. (G07) Their model has two stellar population 
components with different ages, and we adopt here the young component. The age of the old component is 2 Gyr. (L08-1) 
IRAC-undetected LAEs. (L08-2) IRAC-detected LAE s. (P07) The result with a n exponentially decaying star formation history. 
In Figure [T2l we show the results shown in Figure 4 of lCastro Ceron et al. I teOOSl) . (F09) The result with an exponentially 
decaying star formation history. SFRs are calculated from model SEDs. (L07) In Figure 1121 we show their results based on the 
assumption of constant star formation history and 0.005Zq. (C05) They calculate the SFR from the Htv luminosity. 



4-1.1 MlB 



My versus Mis 



Figure [7] plots in large symbols rest-frame UV-to-optical 
colour (A/1500 — My) against UV(= 1500A) absolute mag- 
nitude (M1500) for our LAEs. For conservative discussion, 
rest-frame 1500 A and V-band absolute magnitudes of our 
sample are derived from the best-fit SEDs, not from ob- 
served 7?- and if -band photometry (see details in Section 
133) . 

For the 11 AT-detected LAEs, nine are distributed in a 
similar region in the M1500 — My versus M1500 plane to the 
LBGs. The remaining two are faint in M1500 and extremely 
red (M1500 — Mv > 3). Indeed, these are two of the four 
dusty LAEs, Rl and R2. On the contrary, the /^-undetected 
LAEs, which occupy 96 % of our sample, have bluer M1500 — 
My colours and much fainter M1500 magnitudes than the 
LBGs. The existence of if-bright LAEs suggests that not all 
LAEs at z ~ 3 are low-mass galaxies with blue colours and 
that LAEs are a heterogeneous population with w ide ranges 
of stellar mass, age, and/or dust extinction fe.g..lLai et al.l 
I200&1 : iFinkelstein et al.ll2008l , 120091 : IPentericci et al.ll2009f l 

We also compare our LAEs with those in the literature 
at sim ilar redshifts. The IRAC-detected LAEs bv lLai et al.l 
|200Sl ) have similar M1500 -My and M1500 to the LBGs. On 
the ot her hand, the NIR-un detected LAEs bv lGawiser et all 
120071 ) and lLai et all (|200&t ) have as blue M 1500 -My colours 
and faint M1500 magnitude s as our if -undetecte d LAEs. The 
NIR-undetected LAE by iNilsson etldl (|2007m has redder 
Afi5oo — My than the other NIR-undetected LAEs, probably 
because their if -band image is relatively shallow (Table [TJ . 



4-1.2 Stellar Mass versus My and M1500 

The left panel of Figure [5] shows stellar mass as a function 
of rest-frame IZ-band absolute magnitude for our LAEs cal- 
culated from best-fit SEDs. As expected, the if -detected 
LAEs have comparable stellar masses to the if -detected 
LBGs, spanning over Af star ~ 10 9 - 10 10 ' 5 M Q , while the K- 
undetected LAEs are much less massive, with 10 8 - 1O 8,5 M0. 
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Figure 7. UV-to-optical colour (M1500 — My) versus rest-frame 
UV absolute magnitude (M1500) for LAEs at z = 3.1 and 3.7, 
and LBGs at z ~ 3. The red circles and squares indicate the 
ii"-detected LAEs at z = 3.1 and 3.7, respectively, where those 
with spectroscopic confirmation are shown by bold symbols. La- 
bels Rl - R4 indicate iC-detected LAEs with extremely large 
dust extinction (E(B — V)* ~ 0.65). The blue filled circles and 
squares represent the i<"-undetected LAEs at 2 = 3.1 and 3.7, 
respectively. The dotted and dashed lines correspond to the 3<r 
detection limit in the K band for LAEs at z = 3. 1 and z = 3.7, 
respectively. The black squares are LBGs at z ~ 3 jShaplev et al] 
l200ll:|Papovich et alj|200ll : llwata et al.ll2005l) , all of which are K- 
detected and spectroscopically confirmed. The blue upside-down 
pentagon and diamo nd represent co mposite LAEs at z ~ 3.1 
dGawiser et alj l2006fi and 2 ~ 3.15 llNilsson et alj|2007l) . The 
red and blue triangles ind icate stacked I RAC-detected and un- 
detected LAEs at 2 = 3.1 [|Lai et al.ll2008t ). 
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Our results broadly agree with previous studies. This figure 
shows that LAEs have a wide range of stellar mass. 

The dashed lines indicate four mass-to-luminosity ra- 
tios. It is found that our isT-undetected LAEs have slightly 
lower M s tar/iv ratios than the average if-detected LAEs 
and LBGs, suggestive of younger ages. 

The right panel of Figure [8] plots stellar mass against 
rest-frame UV absolute magnitude. All but two of the K- 
detected LAEs are distributed in a region similar to the 
LBGs. The two exceptions are Rl and R2. The very faint 
M1500 magnitudes with respect to high M sta r of Rl and R2 
are mainly due to the extremely large dust extinction (see 
Table 0. The Jf-undetected LAEs are much fainter in UV 
wavelengths and much less massive than the LBGs. 



4-1-3 Star Formation Rate versus A/1500 and L(Lya) 

The left panel of Figure [9] plots the star formation rate 
against the rest-frame UV absolute magnitude. For all ob- 
jects SFRs are calculated from the best-fit SEDs. For the 
case of E(B — V)* = 0, SFR approximately correlates with 
L1500 , and we plot t his correlation with the dotted line using 
(|Madau et al.l(l998l ): 



SFR [M yr _1 ] = 



1/1500 [erg s *] 
8 x 10 27 



(2) 



Reflecting the wide range of E(B — V)*, the _R"-detected 
objects are widely scattered above the dotted line. Some 
objects with very low extinction are located on the dotted 
line, while the four objects with extremely large extinction 
(Rl - R4) have SFRs more than two orders of magnitude 
higher than those expected from the observed M1500 ■ Except 
for Rl - R4, the Tf-detected LAEs have similar ranges of 
SFR and Mi 500 to the LBGs. 

The iC-undetected LAEs have SFRs of 1 - 100 M Q yr" 1 
and their offsets from the dotted line are within an order of 
magnitude because of the modest dust extinction. NB570- 
if-undetected is offset nearly an order of magnitude because 
of its relatively heavy dust extinction (Table 0. 

The right panel of Figure |TJ] shows star formation rate 
versus observed Lyman a luminosity. Based on t he relation 
betw een Ha luminosity and star fo rmation rate (jKennicuttl 
1 19981 ) under case B approximation (|Brocklehursdfl97ll ) . star 
formation rate is related with Lyman a luminosity as: 



SFR [M yr" 



9.1 x 10~ 43 L(Lya)[erg s" 



(3) 



The dotted line shows this relation. Star formation rates 
derived from SED fitting are found to be higher than those 
from equation (j3j). This can be explained by several reasons. 
Lya photons can be attenuated by dust, and scattered in 
the interstellar and intergalactic medium. Lya photons are 
resonantly scattered and heavily attenuated especially under 
a homogeneous ISM. Intrinsically, the Lya emissivity is also 
affected by the age, metallicity, and IMF. In addition, the 
L(Lya) of objects without spectroscopic redshifts tends to 
be underestimated due to the triangle shapes of the narrow- 
band filters. NB570-_ftT-undetected lies well above equation 
((3|, probably due to its relatively large E(B — V)*. 



4-1-4 Stellar Mass, Age, and E(B — V)* versus EW(Lya) 

The top panel of Figure [10] shows stellar mass versus 
EW(Lya). The Jf-detected objects are massive and have 
a wide range of EW(Lya) over 20 - 220 A, but mostly in 
the range < 100 A. On the other hand, the 7f -undetected 
objects are less massive and have consistently high values of 
~ 150 A. Considering that the Jf-undetected objects closely 
represent the whole LAE population, we can conclude from 
this figure that massive LAEs with high EW(Lya) are rare. 

There is a weak anti-correlation between stellar mass 
and EW(Lya), if the four red objects (Rl - R4) are ex- 
cluded. This a nti-correlation is qualit atively consistent with 
the finding by IPentericci et all (|2009l ) in their LAE sample 
from z ~ 3.5 to 6 that more massive objects have in general 
smaller EW(Lya). 

This anti-correlation is not due to the Lyman a emis- 
sion of more massive objects being more heavily absorbed by 
dust. For example, NB570-W-55371 has a small extinction 
(E(B - V)* = 0.04) and a small EW(Lya) (20 A), while 
NB503-S-94275 with a moderate extinction (E(B - V), = 
0.11) has a large EW(Lya) of 150 A and the 2 = 3.7 K- 
undetected objects with E(B — V)± = 0.19 have a large 
EW(Lya) of 130 A. Indeed, the bottom panel of Figure 
1101 shows no significant correlation between E(B — V)* and 
EW(Lya). A cause of the anti-correlation may be that the 
geometry of dust in LAEs changes with stellar mass in the 
sense that more massive gal axies have a le ss clumpy distri- 
bution for some reason. As iNeufeldl l|l99lf ) points out, the 
fraction of Lya photons escaping from a galaxy is higher 
for a clumpy cloud distribution than for a uniform distribu- 
tion. Th is effect can explain th e observed SEDs of LAEs at 
z ~ 4.5 (|Finkelstein et alj|2009h . 

The middle panel of Figure \TU\ shows stellar age ver- 
sus EW(Lya). There seems to be a weak anti-correlation 
between stellar age versus EW(Lya). 

In the bottom panel of Figure 1101 the small squares 
indicat e the four averaged va lues of 38 LAEs at z ~ 3.5 
- 6 by IPentericci et all ( 20091 ) . A probable reason for their 
very small E(B — V)* is that they have been selected using 
the LBG technique and thus have blue colours. The crosses 
represent the data for local star forming galaxies given in 
iGiavalisco et all (|l99rj ). In comparison to our LAEs, these 
local galaxies have on average smaller EWs. This might sug- 
gest a more un iform cloud distribution for local galaxies, if 
INeufeldl (|l99ll ) picture is correct. 



4-1.5 Lya escape fraction versus E(B — U) gas 
We define the Lya escape fraction / osc (Lya) by: 



/esc (Lya) 



L ob5 (Lya) [erg s 
Lint (Lya) [erg s~ 



(4) 



where L int (Lya) [erg s" 1 ] = 1.1 x 10 42 SFR [Mqvt" 1 ], equiv- 
alent to equation @, is the intrinsic Lya luminosity com- 
puted from the SFR on the assumption of case B. Since 
the SFR is derived from our SED fitting, it has been cor- 
rected for dust extinction. Then / esc (Lya) is plotted against 
E(B - U) gas in Figure ITT1 As for our results, E(B - U) gas = 
E(B-V)*/QM (|Calzetti et al.ll2000t ). A strong negative cor- 
relation is seen between the two quantities for our objects. 
The if-undetected LAEs, whose dust extinction is modest, 
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Figure 8. Left: Stellar mass versus rest-frame V-band absolute magnitude. All symbols are the same as in Figure [7] Black dashed 
lines correspond to four mass-to-luminosity ratios normalised by the solar value. Right: Stellar mass versus rest-frame UV absolute 
magnitude. All symbols are the same as in Figure [7] 
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Figure 9. Left: Star formation rate derived from SED fitting, plotted against the UV absolute magnitud e. All symbols are the same 
as in Figure [7] The dotted line shows the relation between the UV luminosity and the star formation rate llMadau et al.|[T998h , Right: 
Star formation rate derived from SED fitting, plotted against the Lya luminosity. All symbols are the same as in Figure [jj The dotted 
line shows the relation between the Lyman a luminosity and the star formation rate. 



have relatively large / csc (Lya) of ~ 0.1 - 1, while the K- 
detected LAEs, whose E(B — V) sas can be as large as 1.5, 
have / GSC (Lya) ~ 3 x 10~ 4 - 1. Rl - R4 have the smallest 
escape fractions. 

We discuss implications of this negative correlation. 
Assuming /„ = constant over 1200 - 1500 A, we obtain 
^obs(Lya) = EW(Lya)(c/Ao)I/i50o,obs as: 



L obs (Lya) = EW(Lya)^L 1500 10-°- 4ii(Ao)E(s - ,/) -, (5) 



where A = 1216 A, and R(\ ) = 0.44fe'(A ) ^ 5.27 is 
the selective extinction at 1216 A for Calzetti extinction 
law. Combining this equation with equations ([2]) and (|3]), 
we rewrite equation (|4]) as: 



/esc (Lya) 



1 . x 10 42 Ag 



-2-EW(Lya)-10 



-0AR(X o )E(B-V) E 



3 (6) 



The dashed line in Figure [TTJ corresponds to this equation 
for EW(Lya) = 68A, the value of EW(Lycv) expected for 
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Figure 10. Top: Stellar mass derived from SED fitting versus 
rest-frame Lya equivalent width. All symbols are the same as 
in Figure [jj Middle: Stellar age derived from SED fitting ver- 
sus rest-frame Lya equivalent width. All symbols are the same 
as in Figure [7J Bottom: Dust reddening E(B — V)* derived 
from SED fitting versus rest-frame Lya equivalent width . The 
black small squares are taken from [Pentericci et al ] l|2009h . The 
crosses are star-formin g galaxies in the local universe taken from 
Giavalisco ct al ] l|l996h . The remaining symbols are the same as 
in Figure [jj 
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Figure 11. Lyo escape fraction versus E(B — V) ga s. As for 
our r esults, E(B - V) gas = E(B - V)*/0.44 jCalzetti et al.l 
120001) . The triangles ar e star-forming galaxies in the local uni- 
verse llAtek et al]|2008h. and the black small squares are z ~ 6 
LBGs |VeTb*a^mrieet**all I2008T ) . The remaining symbols are the 
same as in Figure [jj The dashed line corresponds to equation 



is interesting that almost all of our LAEs are on the dashed 
line. Similar to the bottom panel of Figure 1101 which shows 
no significant correlation between E(B — V)* and EW(Lya), 
this results also suggests that in heavily obscured LAEs the 
absorption of Lya photons is not as strong as expected for 
uniform ISM. 

Our result is consistent with the result of lValls-Gabaudl 

who has studied nearby starburst galaxies and 
showed that most of them present normal Lya emission con- 
sistent with case B recombination. This implies that the ef- 
fect of resonant scattering i s not so strong. 

iFinkelstein et al. have studied z ~ 4.5 LAEs con- 

sidering their dust geometries by introducing the q param- 
eter as the ratio of the optical depth for Lya to that for 
UV continuum. Our results suggest that Lya photons and 
UV continuum photons near Lya wavelength are equally at- 
tenuated by dust . This corresponds to q = 1. Interestingly, 
IFinkelstein et all (|2009l ) have found that, although three ob- 
jects have q > 1 and two have q = 1, the remaining nine have 
q < 1. This might mean that the dust effect on Lya photons 
is more important at higher redshift. 

The black sm all squares represent th e results for 11 
LBGs at z ~ 3 bv lVerhamme et all l|200ct l. Although they 
also show a negative correlation, its slope is steeper than 
the dashed line. This means that for a given E(B — V) ga , s , 
Lya photons are more strongly absorbed in comparison with 
LAEs. Similarly, local star-forming galaxies (open triangles, 
lAtek et al.ll2008l ) also tend to have lower / esc (Lya). 



/e S c (Lya) = 1 when E(B — V") gas = 0. In other words, we 
assume here that the absorption fractions of Lya photons 
and continuum photons near Lya wavelength are identical 
irrespective of E(B — V) gas . Clearly, this assumption is not 
valid for uniform ISM, since Lya photons are more heavily 
absorbed before escaping from the galaxy. In this sense, it 



4-1-6 Specific Star Formation Rate versus Stellar Mass: 
Comparison with other high-z galaxy populations 

We plot in Figure [12] the specific star formation rate (sSFR 
= SFR/M s tar) against the stellar mass for our LAEs and 
other high-redshift galaxies selected by various photometric 
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selection methods, in order to compare the mass scale and 
star-formation activity between the high-z galaxy popula- 
tions. Also plotted are LAEs over a wide range of redshift 
taken from the literature. The data plotted are summarised 
below: 

LAEs 

Large blue filled circle and large blue filled square: K- 
undetected (stacked) LAEs in our sample at z = 3.1 and 
3.7, respectively. Large red open circles and large red open 
squares: Jf-detected LAEs in our sample at z = 3.1 and 
3.7, respectively, where those with a spectroscopic redshift 
are shown by bold symbols. B lue filled inverse pen tagon: 
18 stacked LAEs at z = 3.1 (|Gawiser et all l2006h. Blue 
filled diamond: 23 stacked LAEs at z = 3.15 |Nilsson et al.l 
120071). Blue filled pe ntagon: 52 stacked LAEs at z = 3.1 
( Gawiser et al.ll2007h . Blue filled triangle and red filled tri- 
angle: 76 IRAC-undetecte d and 18 IRAC -detected LAEs 
at z = 3.1, respectively |Lai et al.l I2008T ). Red open tri- 
angle s: three IRAC-detected LAEs at z = 5.7 ()Lai et al.l 
l2007ft . Blue open diamonds and red open diamonds: IRAC- 
undete cted and IRAC-detecte d LAEs at z — 4.5, respec- 
tively l|Finkelstein et al1l2009l). Blu e open pentagons: nine 



LAEs at z ~ 5 (IPirzkal ct al. 2007). Red open hexagon: a 



LAE at z = 6.56 (|Charv et al.ll2005h 



LBGs 

Small open squares: /("-detected LBGs with spectroscopic 
redshifts at z ~ 3 (Shapley et al. 2001: K s < 24.3, Pa- 
povich et al. 2001: K s < 23.8, Iwata et al. 2005: K < 25.4). 
The distribution of their stellar masses has a lower cutoff 
at ~ 10 9 M Q , which roughly corresponds to the limiting K 
magnitudes. 

DRGs 

Small op en circles: DRGs at z ~ 2-3 selected by J — 
K > 1.3 (jvan Dokkum et al.ll2004l ). which are thought to be 
either passively evolving massive galaxies or dusty starburst 
galaxies. Since they are all /("-bright, their stellar masses are 
higher than ~ 10 n M Q . 

SMGs 

Asterisks: submillimetre galaxies (SM Gs) at z ~ 2 
l)Borvs et aT1l2005l : IChapman et al.ll2005h . which are thought 
to be massive, dusty starburst galaxies. 

Galaxies with Photometric Redshifts 

Black and magenta dots: /-band and Jf-band selected 
galaxies at 2.5 < z < 3.5, respectively, whose redshifts 
are estimated from the photometric redshift technique 
(|Feulner et alj|2005l ^l . 

We first focus on the stacked LAEs (blue filled symbols), 
which should closely represent the whole LAE population. 
The stellar masses of the stacked LAEs spa n over 10 8 - 
10 9 M p) (the individually measured LAEs in IPirzkal et al.l 
(|2007h have slightly lower masses), while those of the LBGs, 
DRGs, SMGs, and photo-z galaxies are 10 9 - 10 n M Q , 10 11 
- lO 12 Af , 10 11 - 10 12 M Q , and 10 8 5 - 10 n M Q , respectively. 
Thus, LAEs are the least massive galaxy population among 
those discussed here. This implies that a significant fraction 



6 http : //www.pik-potsdam. de/members/f eulner/ 
research/cosmology/electronic-data 



of low- mass galaxies beyond the limits of broad-band surveys 
have Lyman a emission strong enough to be detected in 
narrow-band surveys. Note, however, that this does not rule 
out the existence of low-mass galaxies with weak or without 
Lyman a emission. 

The sSFRs of the stacked LAEs span over two orders of 
magnitude, from 10 -9 to 10~ 7 yr _1 , although each measure- 
ment has large uncertainties. This range of sSFR is similar 
to that of the LBGs, implying that these two galaxy popu- 
lations have similar growth rates of stellar mass in spite of 
the large difference in M s t a r- In contrast, the stacked LAEs 
appear to have higher sSFRs than the DRGs, SMGs, and 
photo-z selected galaxies. Since the shaded region has not 
been explored, the existence of low-mass galaxies with very 
low sSFRs is not ruled out. 

We then discuss the K- and IRAC-detected LAEs. Our 
/^-detected LAEs are distributed in a similar region to the 
LBGs (all are /("-detected), except for the f our massive LAEs 
with v ery high sSFRs labelled as Rl - R4. iFinkelstein et all 
(2009) 's IRAC-detected LAEs are distributed in a simi- 
lar region to the photo-z galaxies, with relatively low stel- 
lar masses, whi ch reflects the d eep IRAC Chi limit. On 
the other hand, lLai et al.l \200H ) IRAC-detected LAEs at 
z = 5.7 are as massive as the LBGs due to the relatively shal- 
low limiting magnitude. Overall, the K- and IRAC-detected 
LAEs are not discriminated from the LBGs and photo-z 
galaxies in the sSFR versus M sta r plane. 



4.2 Properties of Red LAEs 

As seen in the previous subsection, the four red LAEs (Rl 
- R4) behave quite differently from the other LAEs in most 
of the two-parameter spaces. They have large stellar masses 
and high SFRs for relatively faint UV absolute magnitudes 
(Figure^. They deviate from a negative correlation between 
M s tar and EW(Lya) seen in the majority of LAEs (top panel 
of Figure 1 10 p . They have large EWs in spite of extremely 
large E(B — V)* (bottom panel of Figure [To"J) . Finally, they 
have much higher sSFRs than other galaxies with similar 
stellar masses (Figure [12)). These results collectively suggest 
that they are not just massive LAEs, but may belong to 
a special class of galaxies involving obscured starburst. In 
this subsection, we further investigate their properties in 
conjunction with MIPS 24/im data and the IMACS spectra. 



4-2.1 Star Formation Rates based on MIPS Photometry 

All four objects have been imaged with MIPS at 24pm in 
the SpUDS. We find that R4 (NB570-W-59558) is detected 
at ~ 5<t level (rri24 M m = 19.3) while the remaining three are 
fainter than the 3a limit (m24 Mm = 19.8). Note that none 
of the four has a counterpart eit her in the XMM-Newton 
X-ray catalogue (|Ueda et al.ll2008T ) with a detection limit of 
/(0.5 - 2 keV) = 6 x 10~ 16 erg cm " 2 s _1 , or in the VLA 
radio catalogue (|Simpson et al.l 2006T I with a detection limit 
of /(1.4 GHz) = 100/iJy. All four are out of the SCUBA 
Half- Degree Extragalactic Survey (jSHADES: Coppin et al.l 
I2006T ). 

We cannot completely exclude the possibility that the 
MIPS photometry of R4 suffers from confusion, because 
there are two objects near R4 in the optical wavebands with 
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Figure 12. Stellar mass versus specific star formation rate (=SFR/stellar mass). The dashed line s represent SFR = 0.01 - 10 Mq 
yr — 1 from left to right . The blue filled pentagon represents a stacked LAE o f 52 objects at z = 3.1 l lGawiser et al.ll2007l) . The red and 
blue open diamonds are 14 LAEs z ~ 4.5 with and without IRAC-detection jF inkclstcin ct aL 200^, private communication). The blue 
open penta gons are LAEs a t z ~ 5 iPirzkal et al 1 l200l l Castro Ceron et al .1120081 ). The red open triang les are three IRAC-detected LAEs 
at z = 5.7 jLai et al.ll2007t). The red open he xagon is a lensed LAE at z = 6. 56 llCh ary et aL]|20p5 [). Th e star marks are SMGs from 
iBorvs et al, U2005I) and lChapman et al. I J2005h . The small circles are DRGs from Ivan Dokkum et al.l j2004h . The black and magenta dots 
are /-band and ff-band detected photo- z samples respectively jFeulner et alj|2005h . The remaining symbols are the same as in Figure [7] 
There is no object in the region of the high stellar mass and sSFR (shaded region in the top-right corner). The lower-left shaded region 
corresponds to quiescent low-mass objects, almost all of which are missed so far. 



a separation of about 1.8 and 3.8 arcsec, respectively. In the 
following, we assume that possible source confusion does not 
significantly affect the MIPS photometry. 

We use the t empla te S EDs of starburst g alaxies given 
in lLagache et al.l (|2003lf l and lVega et~afl (|2008l rlto estimate 
SFRs from the observed flux density. As will be shown 

below, the SFRs thus estimated have considerably large 
uncertainties because corresponds to near the blue 

end of dust emission f or our objects. We u se one template 
(Lih{Lq] — 10 12 ) from lLagache et all (^OOSt ). and templates 



of 14 starburst-dominated ULIRGiPl from lVega etahl (|200ST ). 
The range of their log (L m [L Q ]) is from 11.39 - 12.24. We 
redshift them at z = 3.1 and 3.7 and find scaling relations 
between ni24 (lm and Lir, for individual templates. We then 
derive SFRs from indivi dual Lir us i ng th e relation between 
SFRir and Lir given in iKennicuttJ (|l998l ): 

SFRm [Moyr -1 ] = 1.7 x 1(T 10 Lir[L q ]. (7) 

We obtain SFRir ~ 4.2 x 10 3 - 2.0 x 10 4 [Moyr" 1 ] (median 
SFRir ~ 7.4 x 10 3 [M yr _1 ] ) for R4 and upper limits for 



7 http : //www. ias .u-psud. f r/irgalaxies/model .php#SED 

8 http : //adlibitum. oat . ts . astro . it/silva/grasil/ 
modlib/modlib . html 



9 IR 12112+0305, UGC 9913, IR 10565+2448, IZW 107, IR 
10173+0828, Arp 299, UGC 4881, IC 1623, UGC 8387, UGC 2369, 
IIIZW 35, IC 5298, Arp 148, UGC 6436. 
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the remaining three of l.lxlO 3 -4.3xl0 3 [M©yr _1 ] (median 
~ 1.5 x 10 3 [A^yr" 1 ]) for z = 3.1, and 2.8 x 10 3 - 1.3 x 10 4 
[Mayr -1 ] (median ~ 4.9 x 10 3 [M Q yr -1 ]) for z = 3.7. The 
variation of SFRir among the 15 templates is thus a factor 
~ 4. 

In the left panel of Figure 1131 we compare the SFR 
derived from the SED fitting (SFRsEDflt) with the sum 
of the SFRs derived from the observed UV and IR lumi- 
nosities (SFRuv+ir)- We do not find inconsistency between 
SFRsEDflt and SFRuv+ir, although the SFRuv+ir of three 
objects is an upper limit and thus the constraint on the 
consistency is not so strong. R4 is detected at 24/im and 
its SFRuv+ir is comparable to those of SMGs and DRGs. 
These results support our inference that the four red LAEs 
are massive, dusty starburst galaxies. They can be detected 
with SCUBA2/JCMT and Herschel. 

In the right panel of Figure 1131 the ratio of SFRir 
to SFRsEDflt is shown. The dotted line corresponds to 
SFR tota i = SFRir, i.e., all of the intrinsic UV emission is 
re-emitted from dust. Only R4 is detected at 24pm and has 
SFRm/SFRtotai — 0.2, while the other three have upper lim- 
its of 0.5 - 3. This result may suggest that the contribution 
of the obscured star formation to the total star formation is 
not necessarily very large. 

We make median-stacked MIPS 24pm images of the 
LAEs fainter than K(3a) at each redshift (211 for z = 3.1 
and 70 for z = 3.7). We find that both are fainter than 
m24/im(3cr), obtaining upper limits of SFRir = 1.0 x 10 2 
(5.8 x 10 2 ) [M yr _1 0for the z = 3.1 (3.7) LAEs. The re- 
sults are plotted in blue symbols in Figure 1131 We cannot 
place strong constraints either on the consistency between 
SFRsEDflt and SFRuv+ir (left panel) or on the fraction of 
obscured star formation (right panel). 

Note that we cannot completely exclude AGN- 
contaminated objects from our sample. In particular, the 
MIPS detection for R4 may imply that it has a significant 
amount of AGN contribution. In this case, the SFRs calcu- 
lated above are overestimated. 



4-2.2 Velocity Widths of Lyman a Emission 

As described in Section 12.51 we performed spectroscopy of 
two red LAEs (Rl and R2) with IMACS on Magellan, 
and confirmed them to be real LAEs. We obtain a Lya 



line width of hfwhm = 205 ±129 km 



for Rl, and 



wfwhm = 629 ± 201 km s~ for R2, after correction for 
instrumental broadening on t he assumption t h at th eir in- 
trinsic profiles are Gaussian. IVenemans et al. (2005) have 
reported that the ty pical FWHM of the Lya line of LAEs is 
340 km s -1 (see also Taniguchi ct al. 2005; Kashikaw a et al.l 
l2006r , iTapken et al.ll2007t ). The velocity width of Rl is sim- 
ilar to this value. On the other hand, R2 has twice as large 
a value, and i t might suggest that this object involves gas 
outflows (e.g., iTaniguchi fc Shiovall2000t ) . that it is a merger 
of two objects with a large relative velocity, or that it has a 
heavily obscured AGN. 



10 These are median values of the cases with the 15 SB- 
dominated ULIRGs. The range of the upper limits of SFRir are 
73 - 2.9 x 10 2 [A/oyr^ 1 ] for z = 3.1, and 3.3 X 10 2 - 1.6 x 10 3 
[Moyr" 1 ] for z = 3.7. 



4-2.3 Comparison with Other Red LAEs and Local 
Starburst Galaxies 

We compare in Figure [TJ] the SEDs of Rl and R2 (both 
have a spectroscopic redshift) with those of the thr ee lo- 
cal starburst galaxies Arp 220, M82 (|Silva et alJll99ct ). and 
IR10565+2 448, which has n o AGN contribu tion to its IR 
lumin osity (|Vega et all 120081 1. and Mrk 231 (jPolletta et all 
l2007f ). which is considered to be a local starburst galaxy with 
a heavily obscured AGN. Also plotted are tw o SMGs with 
strong Ly a emission recently discovered by ICapak et al.l 
IMPOST ) and lCoppin et all (|2009h ■ 

Over the rest-frame UV to NIR wavelengths, Rl re- 
sembles Arp 220 and IR10565+2448 more closely while R2 
resembles Mrk 231, and both are almost one order of mag- 
nitude brighter than M82. This might indicate that Rl is a 
pure star-forming system, while R2 is a star-forming galaxy 
with an obscured AGN. The SEDs of Rl and R2 including 
upper limits at the MIPS and VLA wavelengths are con- 
sistent with them being dusty starburst galaxies. The two 
SMGs are even brighter than our objects, Arp 220, Mrk 231, 
and IR10565+2448, and have significantly bluer SEDs. 

Table [TO] summarises their stellar populations. Rl and 
R2 are found to have similar stellar masses to Arp 220, 
IR10565+2448, and the two SMGs (For the z = 4.55 ob- 
ject we adopt the values for the diffuse component); M82 
is ~ 100 times less massive. The E(B — V)* values of Rl 
and R2 are close to that of Arp 220, and the other objects 
are less obscured. Ages are not strongly constrained, but all 
objects appear to be young. The SFRs of Rl and R2 are 
several times higher than that of Arp 220 and are compara- 
ble to that of the z = 4.55 SMG. Thus Table IT01 leads us to 
concluding that Rl and possibly R2 are high- z counterparts 
of present-day ULIRGs. 

The number density of the red LAEs (Rl - R4) in our 
sample is ~ 2 x 10~ 6 Mpc~ 3 at z = 3.1, and ~ 7 x 10~ 6 
Mpc -3 at z — 3.7, respectively. These values are about one 
or der of magnitude low er than that of SMGs at z ~ 2 - 
3 j Chapman et al.ll2005h . If red LAEs and SMGs are both 
massive dusty starburst galaxies, then the difference in the 
number density would imply that dusty starburst galaxies 
emit Lya photons in only a short time in their starburst 
phase. In that time the distribution of dust in the galaxy 
might be clumpy so that Lya photons can escape relatively 
easily. 

Note that we cannot rule out the existence of obs cured 
AGNs in our red LAEs. Actually, iGeach et all (|2009h have 
reported that all AGN in the Lya blobs in the z = 3.09 
proto-cluster in the SSA22 field appear to be heavily ob- 
scured. If the contribution from AGNs to the SEDs is signif- 
icant, our discussion about their stellar populations will have 
to be modified. Even in such cases, however, our finding of 
Rl - R4 demonstrates that Lya surveys can detect galaxies 
with very red UV-to-optical SEDs, which are difficult to be 
identified by other photometric selection methods. 



4.3 Contribution of LAEs to the Stellar Mass 

Density and the Cosmic Star Formation Rate 
Density 

We estimate the contribution from LAEs (both _R"-detected 
and J\-undetected) to the cosmic stellar mass density at 
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Figure 14. Multiwavclcngth SEDs of two spectroscopically confirmed red LAEs (black squares): NB 503-N-21105 (Rl; left) and NB570- 
N-32295 (R2; right ). Also plot t ed are the SEDs of Arp 220 (solid line), M82 (dotted li ne), both from lSilva et al.l jl998ft . I R10565+2448 
(dash ed line) from Ivega et al.l ||2008| '). M rk 231 (chain line) from iPolletta et al.l J2007h . Lya-emitting SMGs at z = 4.87 l lCoppin et al.l 
|2009| . magenta circles) and at z = 4.547 l lCapak et al,|[2008l . green triangles). The inner panels are a zoom up over the wavelength range 
from UV to near-infrared. 



2 = 3.1 and 3.7, by dividing the total stellar mass of LAEs 
by the comoving volume searched by each narrow band. The 
total stellar mass at each redshift is defined as the stellar 
mass of the stacked LAEs multiplied by their number, plus 
the sum of the JsT-detected LAEs' masses. We obtain 1.4 x 
lO 5 [M Mpc~ 3 ] for z = 3.1 LAEs, and 5.2 x 10 5 [Af Mpc~ 3 ] 



for z — 3.7 LAEs. These values can be regarded as lower 
limits of the contribution to the cosmic stellar mass density. 

iGrazian et al.l (|2007h have obtained the stellar mass 
densities of LBGs and DRGs at z ~ 3 by integrating their 
stellar mass functions, to be ~ 1 x 10 7 [MqMpc -3 ] for LBGs, 
and ~7x 10 6 [M o Mpc~ 3 ] for DRGs. The stellar mass den- 
sities of our LAEs are 1 - 3 % of that of LBGs+DRGs. 
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Table 10. Comparison of stellar populations of dusty star burst galaxies 





M * 


/ . ! / J — V jir 




or [ i 


t\ .cicrcnc c 


_t\Grn£LFK 




[Mq] 


[mag] 


[Myr] 


[Mq yr" 1 ] 






NB503-N-21105 (Rl) 


2.7 x 10 10 


0.70 


4.8 


5.6 x 10 3 


This Study 




NB570-N-32295 (R2) 


4.0 x 10 10 


0.67 


25 


1.7 x 10 3 


This Study 




Arp 220 


2.5 x 10 10 


0.6 


(50) 


5.8 x 10 2 


Silva et al. ( 1998) 


1 


M 82 


2 x 10 8 


0.28 


(50) 
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Silva et al. ( 1998) 


1 


IR10565+2448 


2.3 x 10 10 
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Vega et al. C2008) 


2 
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0.37 


40 


5 x 10 2 


Coppin et al. (2009) 


3 


z = 4.55 SMG diffuse 


2 x 10 10 


0.37 
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4.2 x 10 3 


Capak et al. (2008) 
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1 They have used a library generated with GRASIL jSilva ct al,lll998T ) to make model SEDs. The stellar 
mass shown here is the gas mass which was converted into stars during the most recent burst (see below). 
We calculate the colour excess E(B — V)* from the optical depth in the B band they have reported 

(t b ~ 2.8 for Arp 220, and t b ~ 1.3 for M82), using ICalzetti et al.l (§000) extinction curve. For both 
objects, they assumed that the SFR increased in the first 3 Gyr and then smoothly decreased until an age 
of 13 Gyr, followed by a starburst which occurred at 5 X 10 7 years ago. The SFRs shown here are the ones 
averaged over the last 5 X 10 7 years. 

2 They have used a library generated with GRASIL llSilva et al.lll998l) to make model SEDs. The assumed 
an exponentially decaying starburst. For age and stellar mass, we adopt the age of the burst (age]-,) and the 
mass of stars created during the burst (M^). The SFR we adopt is the one averaged over the period of the 
burst. 

3 They have used the stellar population model of lBruzual fc Ch ariot (199j|), the lCalzetti et aL I ll2000l ) 
extinction law, and either a constant or single- burst star formation history, assuming Z = Zq. They have 
derived its stellar mass from an uppe r limit of the rest-f rame i\-band luminosity and the 
mass-to-luminosity ratio obtained bv lBorvs et al.l ||2005| ). They have reported t he dust extinct i on am ount 
in the V band to be Ay ~ 1.5. We calculate its colour excess E(B — V)* using ICalzetti et al. I ll2000l ) law. 
They have estimated the SFR by dividing the stellar mass by an age of ~ 100 Myr. 

4 Since this object has a UV bright knot and a diffuse component, they have analysed separately the knot 
and the diffuse com ponent. We adopt here the values for the di ffuse c omponent. They hav e used the stellar 
population model of lBruzual fc Charlotl ||2003|) . ISaTpeterl dl955l ) IMF. ICalzetti et al.l ifeoOOl) extinction law, 
and a single-burst star formation history, assuming Z = Zq . They have reported t he dust extinction in the 
V band, and we calculate the colour excess E(B - V)* using ICalzetti et al.l l|2000l) law. 



Similarly, we estimate the contribution from LAEs 
(both 7^-detected and JsT-undetected) to the cosmic star for- 
mation rate density, where the total star formation rate 
is computed as the star formation rate of the stacked 
objects multiplied by their number, plus the sum of the 
Jf-detected LAEs' star formation rates. We obtain 8.2 X 
10" 3 [A'/ Q yr _1 Mpc- 3 ] for z = 3.1 LAEs, and 1.3 x 10 _1 [M Q 
yr~ 1 Mpc -3 ] for z = 3.7 LAEs. The cosmic star formation 
density due to all galaxies is estimated to be logp, ~ 0.18 
[M© yr _1 Mpc" 3 ] at z = 3.1, and log p , ~ 0.14 [M Q 
yr^Mpc -3 ] at z = 3.7 (jHopkins et al]|2006l ). Thus the lower 
limit to the contribution from LAEs is about 0.04 - 1 %. 

It is interesting to note that for our LAEs both the 
cosmic star formation rate and the stellar mass density are 
higher at z — 3.7. This may mean that at higher redshift a 
higher fraction of galaxies are actively forming stars, or that 
the fraction of galaxies with Lyman a EWs large enough to 
be selected as LAEs increases with redshift. 



5 SUMMARY AND CONCLUSIONS 

In this paper, we have investigated the stellar populations 
of LAEs at z = 3.1 and 3.7 found in 0.65 deg 2 of the 
SXDF, based on deep, rest-frame UV-to-optical photome- 
try from three surveys made with Suprime-Cam/Subaru, 
WFCAM/UKIRT, and IRAC/Spitzer. Among a total of 302 
LAEs (224 for z = 3.1 and 78 for 2 = 3.7), only 11 (or 4%) 



are detected in the K band, i.e., brighter than K(3a) = 24.1 
mag (AB). Among the 11 /("-detected LAEs, eight are spec- 
troscopically confirmed. In our stellar population analysis, 
we treat if-detected objects individually, while we stack K- 
undetected objects at each redshift to derive an average 
SED. Since the vast majority (96%) are undetected in K, the 
SEDs constructed from stacking should closely represent the 
whole LAE population at each redshift. We treat the LAEs 
at z = 3.1 and 3.7 collectively as objects at z ~ 3, without 
discussing possible evolution between the two redshifts. Our 
LAE sample, based on deep optical and near-infrared data 
over a wide sky area, enables us not only to place strong 
constraints on the average properties of LAEs by stacking 
of many objects, but also to study rare, massive LAEs visible 
in K. We fit stellar population synthesis models to the multi- 
band photometry of T^-detected and K- undetected (stacked) 
objects, to derive their stellar masses, star formation rates, 
ages, and dust extinctions. We assume a constant star for- 
mation rate and fix the metallicity to Z = 0.2Zq. 
Our main results are as follows: 

(i) The if-detected objects have stellar masses of M sta r ~ 
10 9 - 1O 1O ' 5 M , while the A"- undetected objects are 1 - 2 or- 
ders of magnitude less massive, with M stal - ~ 10 s - 1O 8 5 M . 

(ii) The /("-detected objects are distributed in a similar 
region in the M sta r versus My plane to spectroscopically 
confirmed Lyman-break galaxies (LBGs) at z ~ 3, while the 
K- undetected objects have lower M sta r/£v ratios and bluer 
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A^i500 — My colours than the A-detected ones (Figure[7]and 
the left panel of Figure [5} . 

(iii) The star formation rates (SFR) of A-detected ob- 
jects span a wide range of 10 - 10 4 M© yr _1 , while those of 
A'-undetected objects fall between 1 - lOOM© yr _1 . 

(iv) There could be a bimodality in the age distribution 
among the A-detected objects (Figure [6]). The ages of the 
A'-undetected LAEs are around the median value of the ages 
of the A-detected LAEs. 

(v) Four of the A-detected objects have extremely large 
E(B — V) + , and thus their SFRs derived from rest-frame 
1500 A continua and from Lyman a luminosities are both 
underestimated by more than two orders of magnitude 
(Figure [9}. On the other hand, the dust extinction of A'- 
undetected objects is modest. 

(vi) The Lyman a equivalent width (EW) weakly anti- 
correlates with AAtar, except for the four dusty objects hav- 
ing large EWs (Figure [TO")) . No significant correlation is seen 
between EW(Lya) and E(B - V)* (Figure [TO)) . 

(vii) The Lyman a escape fraction decreases with E(B — 
V) gaa . At a fixed E(B — U) gas , the escape fraction of our 
LAEs is higher than those of LBGs and local star-forming 
galaxies (Figure [TTjl . 

(viii) The stellar masses and the specific star formation 
rates (sSFR) of LAEs are compared with those of LBGs, 
DRGs, SMGs, and galaxies with photometric redshifts of 
2 p hot ~ 3. It is found that the LAE population as a whole is 
the least massive among these galaxy populations, but with 
relatively high sSFRs (Figure [12}. 

(ix) Four of our A-detected LAEs with the reddest SEDs, 
two of which are spectroscopically confirmed, are heavily ob- 
scured with E(B — V)i, ~ 0.65, and their continua resemble 
that of Arp 220 and the local ULIRG (Figure[TJJ. They have 
very high sSFRs in spite of their large stellar masses (Fig- 
ure 11 2 p , and one of them are detected in the MIPS 24/im 
image, which supports the SED-fitting result of SFR. These 
facts suggest that the red LAEs are massive, dusty starburst 
galaxies. Our study demonstrates that wide-field Lyman a 
surveys can detect such dusty starburst galaxies. 
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